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M.  Greg  Forest  and  Peter  Mucha,  University  of  North  Carolina  at  Chapel  Hill 

In  this  project,  we  conducted  a  systematic  investigation  on  active  liquid  crystal  flows  and 
flowing  polymer  nano-composites  including  studies  of  nonlinear  phenomenon  in  active 
magnetic  mcrobead  rheology,  detailed  analyses  and  simulations  of  active  liquid  crystal 
models  in  thin  films,  free  surface  geometries,  and  the  channel  geometry,  applications  of 
active  liquid  crystal  models  to  complex  biological  systems,  numerical  algorithm 
development  for  multiphase  fluid  flows,  network  analysis  and  simulations  of 
nanocomposite  systems.  In  the  study  of  active  liquid  crystals,  we  explored  spatial- 
temporal  structures  using  the  two-scale  kinetic  model  by  mapping  out  the  dynamics  in  the 
phase  space  consisting  of  the  active  parameter  and  the  strength  of  active  particle-particle 
interaction.  In  the  meantime,  we  used  a  reduced  order,  continuum  active  polar  liquid 
crystal  model  to  study  the  robustness  of  the  structures  and  their  genesis  in  relation  to  the 
inherent  instability  in  the  active  liquid  crystal  model.  Detailed  studies  are  conducted  with 
respect  to  channel  flows  and  free  surface  filament  flows  which  have  direct  relevance  to 
various  engineering  and  biological  applications.  Network  models  are  brought  in  to 
analyze  nanocomposites  transport  properties  and  network  properties  of  various  complex 
networks. 

New  numerical  algorithms  are  developed  for  a  host  of  multiphase  phase  field  based 
hydrodynamic  models,  in  which  a  new  method  called  Energy  Quadratization  or  EQ  was 
invented  to  systematically  linearize  transport  equations  derived  from  the  generalized 
Onsager  principle.  These  new  methods  are  then  applied  to  efficiently  integrate  the 
governing  system  of  equation  for  the  active  matter  systems  are  continuously  improved. 

We  present  experimental  data  and  numerical  modeling  of  a  nonlinear  phenomenon  in 
active  magnetic  microbead  rheology  that  appears  to  be  common  to  entangled  polymer 
solutions  (EPS).  Dynamic  experiments  in  a  modest  range  of  magnetic  forces  show  (1)  a 
short-lived  high  viscosity  plateau,  followed  by  (2)  a  bead  acceleration  phase  with  a  sharp 
drop  in  apparent  viscosity,  and  (3)  a  terminal  steady  state  that  we  show  resides  on  the 
shear-thinning  slope  of  the  steady-state  flow  curve  from  cone  and  plate  data.  This  latter 
feature  implies  a  new  protocol  to  access  the  nonlinear  steady  state  flow  curve  for 
biological  EPS  available  only  in  microliter-scale  volumes.  We  use  the  moment  closure 
form  of  the  Rolie-Poly  kinetic  model  for  EPS  hydrodynamics,  together  with  a  decoupling 
approximation  that  obviates  the  need  for  a  full  three-dimensional  (3D)  flow  solver,  to 
qualitatively  reproduce  this  dynamic  experimental  sequence.  We  thereby  explain  the 
phenomenon  in  terms  of  entangled  polymer  physics,  and  show  how  the  nonlinear  event 
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(acceleration  and  termination  on  the  shear-thinning  response  curve)  is  tunable  by  the 
interplay  between  molecular-scale  mechanisms  (relaxation  via  reptation  and  chain 
retraction)  and  magnetic  force  controls.  The  experimental  conditions  mimic  movement  of 
cilia  tips,  bacteria,  and  sperm  in  mucus  barriers,  implying  a  physiological  relevance  of  the 
phenomenon  and  compelling  further  quantitative  kinetic-flow  3D  numerical  modeling. 

We  have  developed  a  set  of  kinetic  models  and  simulation  tools  to  study  anisotropic 
nano-particle  dispersions  and  active  particle  dispersions  to  generate  spatial-temporal 
structures  that  the  ensemble  can  exhibit.  Using  the  kinetic  models,  we  have  investigated 
the  spatial-temporal  structure  of  the  model  solutions.  In  the  end,  we  have  obtained  the 
phase  diagram  for  the  solutions  in  the  important  model  parameters.  Using  a  coarse-grain 
model  for  active  polar  liquid  crystals,  we  analytically  and  numerically  explored  the 
spatial-temporal  patterns  in  active  particle  systems  and  studied  the  solution  instability  in 
details  in  channel  flows  and  free  surface  liquid  filaments.  In  particularly,  we  investigated 
the  capillary  instability  of  a  falling  active  liquid  crystal  filament  and  discovered  a 
stunning  set  of  instabilities  due  to  the  activities  of  the  liquid  crystal  molecules  or  active 
nanoparticles.  We  anticipate  that  this  will  lead  to  new  applications  of  active  particle 
dispersions  in  material  processing  and  interpretation  for  capillary  dynamics  in  biological 
systems. 

After  quenching/annealing,  we  adopt  the  statistics  in  the  probability  density  function  to 
analyze  the  materials  properties,  primarily  electric  properties,  using  network  models  and 
assessment  tools  developed  for  community  networks.  The  work  performed  builds  on 
kinetic  models  and  corresponding  solution  algorithms  for  processing  conditions  of  nano- 
particle-endowed  films,  coatings  and  membranes.  Drawing  from  databases  generated  by 
these  codes,  we  are  actively  investigating  network-mediated  properties  that  are  not  well- 
modeled  by  homogenization  theory.  We  have  explored  DC  resistive  electrical  response  of 
connected  nanorod  composites  at  particle  volume  fractions  above  the  percolation 
threshold.  These  calculations  are  performed  in  the  high  aspect  ratio  conductivity  limit, 
wherein  only  currents  along  particulates  are  considered.  As  an  extension  of  the  above 
investigation,  and  as  a  means  towards  understanding  the  general  behaviors  observed  there, 
we  have  also  considered  similar  effects  and  observations  obtained  in  the  classic  “random 
resistor”  setting  of  bond  percolation  on  3D  cubic  lattices.  We  have  also  investigated 
dielectric  properties  at  particle  volume  fractions  below  the  percolation  threshold.  In  these 
calculations,  the  associated  network  is  an  all-to-all  capacitance  network. 

New  and  efficient  numerical  algorithms  for  multiphase  fluid  flows  formulated  using 
phase  field  are  developed  to  decouple  the  momentum  transport  equation  from  the  phase 
transport  equations  making  the  numerical  computational  tool  more  robust  and  modular. 
Some  of  these  numerical  schemes  are  applied  to  study  cellular  aggregate  fusions  and 
liquid  jet  dynamics  of  liquid  crystal  polymers. 

We  also  applied  the  active  liquid  crystal  models  and  the  numerical  technology  to  develop 
a  set  of  new  models  and  simulation  tools  for  cell  motion  in  mitosis  and  migration  on 
patterned  substrate.  These  models  are  being  used  in  collaboration  with  applied  physicists 
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and  material  scientists  to  study  the  interaction  between  the  cell  and  the  substrate  material 
in  order  to  design  novel  cell-sensing  and  guiding  materials. 

5  PhD  students  have  been  partially  supported  by  this  award  during  their  graduate  studies. 
They  have  graduated  and  move  on  to  either  continue  their  postdoc  training  or  working  in 
universities  and  industries. 

The  team  has  been  very  productive  during  the  last  four  years.  Papers  published  and 
submitted  in  the  period  are  listed  below: 
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